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Abstract 

We consider heat transfer processes in an ice hockey hall, during operating 
conditions, with a bottom-up approach based upon on-site measurements. De¬ 
tailed temperature data of both the ice pad and the air above the ice rink are 
used for a heat balance calculation in the steady-state regime, which quantifies 
the impact of each single heat source. We solve the heat equation in the ice slab 
in transient regime, and obtain a general analytical formula for the temperature 
profile. This solution is then applied to the resurfacing process by using our 
measurements as (time-dependent) boundary conditions (b.c.), and compared 
to an analogous numerical computation with good agreement. Our analytical 
formula is given with implicit initial condition and b.c., therefore it can be used 
not only in ice halls, but in a large variety of engineering applications. 

Keywords: transient heat conduction, cooling, theoretical models, analytical 
solutions, ice rinks, energy efficiency, energy saving, eigenfunctions expansions 


1. Introduction 

Ice hockey arenas constitute a challenging setup for studies in building 
physics, due to the diverse physical processes, often concurring to each other, 
which take place in the environment. For instance, ventilation and air condi¬ 
tioning in such large and complex systems are an intriguing field of study. 

Even more importantly, these structures consume a very large amount of 
energy, approaching ^ ISOOMWh per year. In particular, refrigeration is re¬ 
sponsible for nearly 43 percent of the total energy consumption of the ice hall 
[Dlllllll]. Process optimization and energy saving have therefore become key 
concepts in research in the latest years. 

Heat and mass transfer processes take indeed a relevant part in the en¬ 
ergy balance inside an ice hall, as they occur both above and inside the large 
ice/concrete slab forming the ice rink. A comprehensive literature about this 
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topic has been recently produced; thermal processes such as radiation from the 
ceiling and heat transfer between the ice surface and the refrigeration system 
have been thoroughly studied, see e.g. BW- 

In particular, most studies focus on larger scales and rely importantly on 
simulations and/or numerical methods [31 HI [HI HOI E], with a few exceptions 
m- Here we address the topic on a smaller scale, i.e. an ice pad element, using 
a bottom-up approach that comprises both accurate experimental data and the 
derivation of an analytical temperature profile in the ice slab. Moreover, we 
choose the resurfacing stages because they are the most complex and energy 
consuming phases of the operational cycle of an ice hockey hall, due to the large 
amount of heat transferred to the ice pad in a relatively short time. This affects 
importantly the efficiency and energy consumption of the refrigerating system, 
as discussed e.g. in [1]. 

The present article is aimed to finding applicable quantitative knowledge, 
which can be used in the actual efforts towards energy costs reduction. In our 
paper [T] , we investigated the refrigerating system and studied energy consump¬ 
tion and optimization possibilities. Through an analytical calculation, we found 
the optimal configuration for a maximal coefficient of performance (COP). Here 
we focus on a smaller scale, and consider only the heat transfer processes which 
affect the ice pad temperature. 

Though obtained for a specific case-study, our measurements can be applied 
to most of ice hockey halls, for instance in the development of new control 
methods for ice rink cooling systems dam HU IS]. Our findings also help in 
developing methods to collect a portion of the heat generated (see Eq.Q) and 
reuse it in the ice hockey hall. Moreover, the temperature profile formula (541 
is very general and can be applied to a whole class of heat conduction processes. 

The present paper is organized as follows: 

Section addresses the experimental setup and measurements. We describe 
the ice hall, the ice/concrete slab and the different phases of resurfacing. Our 
indoor air data provide a stratification mapping at several heights from the 
ice rink, from 4cm to 8m above the ice, together with temperature and heat 
flux measurements at the ice/concrete interface. An according energy balance 
analysis is also performed and validated with the experimental data, whose 
interpolation curves are integrated and give a satisfactory agreement with the 
theory. 

In Section we examine the heat flux on the ice pad, the relative humidity 
(RH) close to the ice and the temperatures at surface and at the ice/concrete 
interface, under the steady-state conditions preceding the resurfacing. Detailed 
heat balance calculations pertain each contribution (such as the radiation from 
the ceiling and the effect of the lighting system) separately. 

We find the largest contribution to be thermal radiation from the ceiling 
(74%), followed by lighting (14%), with excellent agreement between our calcu¬ 
lations and the temperature and heat flux measurement^ 


^This also constitutes a method for checking the data consistency. 
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Figure 1: A portion of the ice rink floor structure considered. 


In Section we derive an analytical formula for the ice pad temperature 
in transient conditions, and verify its consistency in two limiting cases. In 
Appendix A we then apply it to the resurfacing process, using our measured 
data as initial condition and boundary conditions. Our result is compared with 
a numerical calculation, showing good agreement. 


2. Ice hockey hall and temperature measurements 

Good skating conditions occur when the ice surface is smooth, hard and 
slippery enough. Usually the ice temperature at the surface is between -3°C' 
and -5°C'. In any cases the temperature of the ice should be kept under -1°C to 
maintain ice hardness m- To keep the ice smooth and in optimal conditions 
after the wear due to skating, it is necessary to perform periodical maintenance, 
which is usually done by means of resurfacing machines. These devices first 
shave the ice, then brush it and eventually spread a thin layer of new water on 
the surface. 

At each maintenance cycle, 300 to 800 liters of water are used, corresponding 
to a 0.25-0.5mm thick water layer. The water cooling and freezing processes 
generate a sudden increase in the ice temperature, which lowers the cooling 
efficiency of the refrigerating system, thus increasing the operational costs of the 
ice hall. Accurate analysis of this phenomenon is then important to optimize 
the refrigeration process and to achieve overall energy saving. 

The field study discussed in this paper was carried out in August 2013 at the 
Reebok Arena in Leppavaara, Finland. The ice hall has two identical ice rinks 
of size 1624m^ each. There are no major stands in the hall and the ventilation 
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Figure 2: Measurements of temperature and heat flux at the ice/concrete interface. 


air flow rate varies between 3 — 6 /s during the year. The refrigeration system 

is indirect, with ammonium as the refrigerant and ethylene glycol as the brine. 

Consider now a modular component of the ice/concrete slab as in Fig0 ice 
30 mm, concrete slab 150 mm with brine pipes at the middle, EPS-insulation 
100 mm, sand 100 mm with frost pipes and a load bearing concrete slab 250 mm. 
The coefficient of performance for the refrigeration system is ^ 1.6 on the aver¬ 
age [T]. 

We focus on the first resurfacing process in Figj^ corresponding to the peak 
near 12:45. The ice surface has an initial surface temperature Tg ~ —4.5°C, the 
whole ice pad is at —5°C' on the average. During resurfacing, m^ = 450 fcg of 
water at temperature = 40°C' are spread on the ice surface at t = Os. The 
entire process consists of three phases: 

1. cooling of the water layer by contact with the ice pad, convection and 
radiation with the surroundings, 

2. complete freezing of the water by the same processes, 

3. cooling of the new ice from T/reeze = 0°C to ^ —4°C (see Figj^. A 
thermal camera records the surface temperature point wise, taking pictures at 
intervals At = 10 s. Moreover, heat flux and temperature are measured with 
a heat flux plate and a pt-100 temperature sensor installed at the ice/concrete 
interface in the slab, as shown in Figj^ The according data are plotted in FigH 

The physical properties of water are evaluated at the film temperature 7/ = 
20°C'. This implies the specific heat = A.182 kJ/kgK. Since m^ = 4:50 kg 
of water are spread on the ice, the average thickness is only 

a; = = 0.28 mm , (1) 

which is anyway just an indicative value, as the precision of the resurfacing 
machine is not high enough. The size of the ice rink is A = 1624 m^, hfs = 
338 kJ/kg is the water latent heat of freezing and Cp^i = 2.05 kJ/kgK is the 
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Figure 3: Heat flux and temperatures during the water chilling and freezing. 


specific heat of ice at T ^ 0°C. The water chilling and freezing loads are thus 
obtained as follows, 


Qw 

— [?77-^ [/ijs 0) -|- Cp^i 

(0- 

- T,ce)]] 


(2) 

k J 

= 231.21 MJ; q^ = 142.37^ , 




(3) 




Qi 

9i = ^ = 

kJ 

-- 46.4 — , 

Qi 

= rriujCp^wiTw ~ 0) = 75.28 MJ ; 

Q 2 

T 

( 4 ) 

Q 2 

= hf giTiw = 752.1 M J; <12 = 

„ kJ 
= 93.7^ 

kJ 

2.27^, 

( 5 ) 

Q 3 

= Wu,Cpy(0 - Tice) = 3.69 MJ ; 


Q 3 

93 = ^ = 

(6) 


where Tice is the temperature of the frozen resurfacing water under the new 
steady state conditions. The theoretical heat load Eq.([^ is consistent with the 
literature [HIS]. 

Let us verify immediately that the theoretical heat flux in is consistent 
with the measurements at the ice/concrete interface in FigjS The integral 
of the heat flux curve is the total heat transferred during the resurfacing, to 
be compared with Eq.([^. For better precision, we split the curve into three 
contributions: two polynomials plotted in Figj^and Figj^ and a constant value 
qmax = 85.4839 VF/m^, corresponding to the narrow {At ~ 30s) plateau at the 
top of the curve in Figj^ 

If ciA is the heat flux rate for the first contribution and qb the flux for the 
second curve, we obtain the following expressions: 


q^{t) = 3x +0.0106^2 +0.0414t +45.191, (7) 

qB{t) = 3 X - 2 X + 5 X - OMAlt + 85.873. (8) 
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Figure 4: Interpolation curve for Eq. 


Figure 5: Interpolation curve for Eq. 


Integrating the above over the respective time intervals gives 

/•140 ^ j 

qA= qAit)dt = 10.05 ^ , 

Jo 

r2300 j^j 

qB= qB{t)dt = 127.88 — , 

Jo ™ 

to which we add the heat transferred at the peak, namely 


Qmax 


W k J 

= 85.4839 ^ X At [s] = 2.57 ^ . 


(9) 

( 10 ) 


( 11 ) 


Thus the total heat transferred to the ice pad during the three phases of resur¬ 
facing in Eq.(|^ is measured as 

kJ 

fep = 140.49 ^ . (12) 

This is very close to the theoretical mean value in Eq.([^, namely 

k I 

= 142.37 ^ , (13) 

we find indeed 

k 1 

Ag = g™ - qexp = 1.88 —, %(Ag) = 1.32% , (14) 

which is fairly satisfactory and verifies the consistency of the theoretical and 
measured heat loacE 


^This is precise enough for our purposes, notice anyway that one has some freedom in 
setting the upper limit of integration in ( |10[ |. The order of magnitude is anyway what matters. 
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2.1. Indoor air measurements 

Regarding the air temperature and RH, one set of data was obtained at 
0.04m, 0.1m and 0.24m above the ice rink with accuracy ±2% for 0-90% rel¬ 
ative humidity and ±3% for 90-100% RH. The accuracy of the temperature 
measurement is ±0.4°C. Additional measurements were made at 0.6m, 1.2m, 
1.8m and 2.2m above the ice with PT-100 temperature sensors, with accuracy 
±0.5°C. Moreover, the temperature was measured at 5.0m and 8.3m height 
with an accuracy of ±0.5°C. All measurements were done with At = 10s in¬ 
terval and lasted approximately four hours. The probes set at 0.04-2.2m height 
were set up in a metal holder standing above the ice; this was located at one 
corner of the ice rink, Im from the edge of the ice pad. 

A second set of measurements was made in April 2015 to verify and comple¬ 
ment the earlier results. This measurement was done at 0.005m, 0.01m, 0.04m, 
0.1m, 0.24m and 5m height over the ice. The accuracy of relative humidity of 
the probe was the same as in 2013, while for temperature measurements the 
accuracy was ±0.3°C'. The measurements were done at one instance after the 
sensors had became steady. The probes at 0.005-0.24m height were again set 
up in a metal holder over the ice. Two locations of the holder were used, at the 
corner of the ice rink and close to the center of the ice pad, at approximately 
8m from the edge. 

The August measurements for the air temperature at 0.6m-8.3m over the ice 
are plotted in Fig|^ and those at 0.04m-0.24m in Fig0 The air temperature, 
relative humidity and absolute humidity measured in April 2015 are shown in 
Table □ 

The plots show that, as expected, the air temperature above the ice rink 
is strongly stratified, see Figj^ Figj^ and Table The temperature at 0.6m 
height is 2°C and \{)°C at 5m height. Interestingly, by virtue of the mixing 
effect of the ventilation system shooting air at 25°C' at 5m above the ice, the 
highest temperature is at 5m instead of 8.3m (Figj^. Furthermore, the stratified 
air temperature at 0.6-2.2m turns into a more uniform temperature when the 
skaters enter the ice rink at 14:30 and start mixing the air layers. 

The measured relative humidity (RH) and absolute humidity (AH), listed in 
Table suggest that above the ice rinks the RH is uniform at heights 0.005m 
- 0.24m over the ice and distinctly lower farther from the ice pad. The relative 
humidity is higher above the ice than in the other parts of the space, yet it never 
reaches 95-100%. 

The measurements did not imply any significant differences in temperature 
or relative humidity depending on the location over the ice pad (corner or central 
location). 
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Figure 6: Air temperature above the ice during resurfacing and skaters activity. 



Figure 7: Air temperature above the ice rink at different heights, closer to the ice. 


Height over the ice [m] 

T [°C] 

RH [%] 

AH [g/m^] 

5.0 

4.2 

77.2 

5.0 

0.24 

-0.8 

90.2 

4.1 

0.1 

-2.2 

89.6 

3.7 

0.04 

-2.7 

89.3 

3.5 

0.01 

-3.1 

88.1 

3.3 

0.005 

-3.5 

88.2 

3.2 


Table 1: Air temperature, relative humidity (RH) and absolute humidity (AH) as measured 
at different heights from the ice pad surface. 











































































3. Heat flux on the ice track and temperatures before resurfacing 

In this section we use the field measurements to compute some estimates 
of the heat loads on the ice pad before resurfacing. We consider convection, 
condensation and radiation. In this case the heat transfer is steady state in 
good approximation, so this is easy to accomplish. Even though it is difficult 
to obtain very precise values, this shows well the various factors concurring to 
the overall energy balance on the ice hockey rink. 

Consider first Figj^and Figj^before resurfacing, namely before ~12:43. The 
measured heat flux through the ice pad at the ice/concrete interface is on the 
average q = 41.85 IF/m^, and the temperature is T/ = —5.2°C. The ice pad 
thickness \s L = 30mm. We can derive the ice temperature at surface very 
easily, if = ‘2.25W/mK, 


Ts = Ti + ^q = -4.64°C, (15) 

^ice 

which is consistent with the ice surface temperature given in Figj^ 

The heat flux on the ice track surface is the sum of several contributions, 


q — qrad 4“ qconv 4“ qwvcond 4“ qianip 

^rad{.^ceiling 4“ {h(.Qnv 4“ ) 4“ qiamp : 


(16) 


namely thermal radiation from the ceiling, convection and water vapor conden¬ 
sation at the surface, and heat load from the lighting system. For simplicity, 
we neglect the thermal radiation from the vertical walls and from the audience 
stands (they do not give a relevant contribution anyway). 

For convection we use = T(0.04m) = —3.5°C, see Table Accordingly, 
the heat transfer coefficient takes into account both natural convection and a 
correction given by forced convection [20] 

W 

/icon. = 3.41 + 3.55 F = 3.94^-, (17) 

corresponding to F = 0.15 m/s for the air flow right on top of the ice, specifically 
at the height 4cm, as stated above. The convection heat flow is therefore 


Qconv — ^convi^in ^s') — 4.49 


IF 


W 


The radiation heat transfer coefficient is written instead as 

hrad = £12 C’’ (T/gjjj^g 4- T.i^ce.sur) {'^ceiling 4" Tg) = 1.39 ^2]^ 
where TceUmg = 18°C and the resulting emissivity is 

1/1 Ag / 1 


£12 = 


. ^ci \ ^ceiling 


-11 + /: 


Al V £i, 


- 1 


(18) 


(19) 


( 20 ) 
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The view factor is Fd = 0.68, and emissivities Sice.sur = 0.98 and eceUmg = 0.28 
for the ice surface and the ceiling (a load bearing sheet of galvanized steel) are 
respectively used [2]. These give £12 = 0.28. 

One must also take into account the radiative heat transferred to the ice pad 
by the lighting system. The lamps in Leppavaara are metal halide, which implies 
a contribution of 400 W per lamp. The portion of this power that is turned into 
heat is nearly 62% [32]. Using the upper limit for the heat generation for 40 
lamps gives the following contribution: 


Qiamp = 9.92 kW, (21) 

corresponding to the following heat flux, 

= ( 22 ) 

A 

The water vapor condensation heat load is computed via the formula oin 


Qwvcond — ) 


w 




(23) 


where the heat transfer coefficient for condensation hd is calculated from 

Ap 

hd = ITbOhconv , [p] = [atm] 

Ap = (fiinPin - Ps ■ (24) 


Here AT = 
ice surface, 
[T]=[C]) 


Tin — Ts^ and pin = 0.88 is the relative humidity at 4cto from the 
as in Table The saturation pressures are calculated from (here 


Pin = 10® exp ^17.391 — 
Ps = 10® exp I 17.391 - 


6142.83 
273.15+ T, 
6142.83 
273.15 + Ts 


(25) 

(26) 


We thus obtain hd = Q.9W/m?K and qwvcond = 1.03VF/m^. By substituting 
this result into (16), together with Eqs.(18l, (19) and (20), and adding also (22), 
we get 


q — qrad~h qconv qwvcond qiamp — 31.254.49 T1.036.11 — 42.87 


w 


(27) 


that overestimates only slightly the measured value 41.85 The percent¬ 

age of each contribution is listed in Figj^ 
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Figure 8: Contributions to the heat load over the ice pad in steady state conditions, Eq|2^. 


Now we focus on the ice/concrete pad and consider only conduction. The steady- 
state conditions provide for the heat flow inside the slab 


<7 = 


Ts-Tj, 


Ts-Tp 


Rh 


Ri, 


Rr 


(28) 


where Tp is the temperature at the top of the pipes and the thermal resistances 
of the ice and concrete slabs are written as 


Rirp. 


L 


= 1.33 X 10“ 


Km^ 


Rr.nn.c — 


= 1.67 X 10 


-2 


w ’ 

Km? 


W 


(29) 

(30) 


where kice = 2.25 W/mK and kconc = l.SlT/mitT. The temperature at the top 
of the pipes Tp is then estimated from 


(j = 41.85 IC/to^ = 


2 _Ts - Tp 


Rt. 


(31) 


which gives Tp = —5.9°C and can be cross-checked with heat balance inside the 
concrete slab only, 

Ti-Tp 


q = 


Rr 


(32) 


^conc 

This in fact returns q = TV.92Wjm? « 41.85IV/m^. 

Let us finally check the agreement between the thermal camera measure¬ 
ment and the result (151 for the ice surface temperature. The heat transfer via 

qi = — ■ ( 33 ) 


conduction inside the ice pad is written as 

Ts-Tj 


R,. 
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Since this also must be equal to q, namely qi = q = 41.85 Wjrr? because of the 
steady state conditions, substituting and solving with respect to the temperature 
at the ice/concrete interface we get 

Ti=Ts-R^aeq = -^-^°C, (34) 

that is consistent with the measured value indeed. 


4. Analytical temperature profile for the ice pad 


In this section we derive an analytical formula for the temperature profile of 
the ice pad Tice{t,x), which can be used under any specific situation occurring 
in the ice hockey hall. 

The problem consists of solving the heat equation 


du d'^u 

dt dx'^ 


(35) 


where 0 < x < L = 30mm, and aj is the thermal diffusivity of ice, with the 
time-dependent boundary conditions 


u{0,t) = Ts{t), (36) 

uiL,t)=Ti{t), (37) 


and the initial condition 

u{x, t = 0) = fix) = 18.67(0.03 - x) - 5.2, (38) 

which is easily retrieved from the temperature data in the steady state regime. 
It gives indeed Ts(0) = —4.64°C' at the ice surface and T/(0) = —5.2°C' at the 
ice/concrete interface, as in Fig[Tg 

The boundary conditions are in general given by the measurements. In this 
specific case, Ts{t) is computed at the ice surface (at the water/ice interface) 
and T/(t) at the ice/concrete interface. They both are illustrated in Figj^ 

To obtain the analytical form of 7s(t), we interpolate the temperature of 
the ice pad at the surface. The overall trend, including the entire curve from 
the beginning of resurfacing to t = 620s, is clearly logarithmic. It gives the 
equation 

Ts{t) = -0.641 ln(t) -t 0.4016 t[s] < 620 . (39) 

If instead we consider e.g. only the first minute, the surface temperature be¬ 
comes 

rs(t) = 2 X 10-43-0.00151^ - 0.0065t- 1.0633, t[s]< 60. (40) 

These interpolations are given in Figs[^ and The temperature raise in the 
last 10s is due to the third order polynomial (this is non physical, still it remains 
within the experimental error and does not affect the result sensibly). 
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Figure 9: Interpolation of the surface tem¬ 
perature Ts{t) over t < 620s. 


Figure 10: Interpolation of the surface 
temperature Ts{t) for t < 60s. 


The formula for the ice-concrete interface temperature, in the first approx¬ 
imation, is also a simple third-order polynomial obtained from the field mea- 
surements (Figj^, 

Ti{t) = 6 X - 3 X + 4.3 x - 5.2071, (41) 


which is plotted in Figfl^ We will apply these in the Appendix. 


To solve the Cauchy problem given by Eqs.([35|) 
adopt the method of Eigenfunctions Expansions 
impose the following Ansatz 


(|36j), (|37|) and (|38|), we 
113^ as follows. First 


that implies 


w{x,t) = Ts{t) +x 




w{0,t) = Ts{t), 
w{L,t) = Ti{t). 

Now define the temperature profile u(x,t) as the sum 

u{x, t) = w{x, t) + v{x, t), 


(42) 


(43) 

(44) 


(45) 


and substitute this expression in the Cauchy problem. We obtain 


— ^I^XX 1 


(46) 


where we have simplified the derivatives notation, with the Dirichlet boundary 
conditions ?;(0, t) = v{L, t) = 0 and the initial condition v{x, 0) = f{x) — w{x, 0). 
Now use the Eigenfunction expansion 

OO 

v{x,t) = '^Vn{t)s\a{XnX), (47) 

n—1 
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Figure 11: Initial condition f{x). Figure 12: Interpolation of the interface 

temperature Tj(t) for t < 60s. 


to separate space and time dependence. The eigenvalues and eigenfunctions 
associated to the Dirichlet b.c. are 


( TITT \ 

— j , n e N ; Xn(x) = sin (A„x) 
S{x,t) = -wt = - TsW) (I) - ts{t). 


Therefore we obtain a first order linear ODE 


^ + aiX^Vn = Sn(t ), 


(48) 

(49) 

(50) 


where 


rL r 


Suit) = - 


-Ts{t) - -fs{t)j sm(J^^dx 


2 

niT 


f sin? 


\ niT 


- 1 Ts{t) - 


( sin? 


\ TVK 


— cosnTT I Tlii) 


(51) 


with an integrating factor F(t) = Eq.(50l can then be integrated to give 

the following solution, 

x(x,t) = ^|y sin(A„x), (52) 

where the coefficients 


C-n. — 


^ {/(d- (T/(0)-rs(0)) (I) +rs(0) Isin^*^) dx, (53) 

are given by the initial condition u(0,x) = /(x). 

Putting everything together, we can now write the analytical temperature 
profile for the ice pad Eq. (45) as 

u{x, t) = Tice{x, t) = [Ti{t) - Ts{t)] + Ts{t) 

+ ^ I/" sin(A„x), (54) 

n=l Oo J 
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that constitutes a novel result of this paper. This is a completely general for¬ 
mula, with implicit initial condition f{x) and b.c. T/(t) and Ts{t). 

The above clearly reduces to the ice temperature Ts{t) at the surface for 
a: = 0, while it is slightly less immediate to verify that ( [54| is consistent with 
the temperature at the bottom of the ice pad T/(t) for x = L. In this case we 
get 


u(L, t) = T,,e{L, t) = Ti{t) = [Ti{t) - Ts{t)] + Ts{t) 

+ ^ I /" -I-sin(A„L). (55) 

Recall now that 

A„ = ^ , (56) 

which implies 


CXD . / •* j- \ OO 


Sin 71TT 


= 0 , 


(57) 


since each term in the summation is identically zero. Regarding the term con¬ 
taining the integral, 


\ 2 pt 

2 y^ ^ " sin (An£)e~”^'^"* / 

I Jo 


OLiXt, T 


Ts{t) oc E nsinnTre 




= 0, (58) 


again because sinnar = 0, Vn G JV. So we get an identity 0 = 0, as required. 


Fisically, the terms in brackets in Eq.(54) give the transient state correction 
to temperature, that depends on the history of the process (via the t-integral) 


and on the initial and boundary conditions by virtue of Eqs.(51) and (53). This 


is plotted in Eig |A.18|for our specific example calculated in the Appendix, where 
the initial condition (381 and the boundary conditions (40) and (41) pertain the 
resurfacing process. 

In 


The according analytical temperature profile (54) is given in Figure 13 


Fig[^ this is compared to a numerical solution by a Finite Element Model 
(EEM), which is plotted in Figure [l4| Details of the analytical computation are 
discussed in the Appendix. 

In conclusion, Eq.( |54[ ) provides a general formula for the temperature profile 
T(x, t) at any point x of the ice pad, at any generic time t corresponding to any 
situation in the ice hockey hall (closing hours, resurfacing or skaters’ activity). 
The practical problem of studying the complex physical processes over and 
under the ice pad is here avoided by virtue of the boundary conditions (i.e. 
experimental data), which encode all the involved phenomenology in a very 
simple analytical form. 
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Figure 13: Temperature profile in the ice pad, analytical solution. 



Figure 14: FEM solution for t = 0...60s. 



Figure 15: Temperature profile after t = 30s, analytical (solid) versus FEM (dashed). 
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5. Conclusions 


In this paper we have considered thermophysical processes in an ice hockey 
hall during standard operation hours. Detailed heat flux, air and ice temper¬ 
ature and relative humidity data are provided and discussed in a quantitative 
analysis with specific focus on the maintenance (resurfacing) phase. An energy 
balance calculation shows the different contributions to the heat load on the ice 
rink, and we find an analytical formula for the temperature inside a medium 
which is subject to conduction with time-dependent boundary conditions. We 
apply this formula to the case at hand, and obtain a good agreement with an 
analogous numerical computation. 

Outside the occupation hours, we find a strong temperature stratification 
in the air above the ice rink, which is compromised once the skaters start their 
activity. The effect of the ventilation system, delivering air at 25°C at about 5m 
above the ice, is instead independent of the skaters, yet it is critically affecting 
the air layers temperature. As it is shown in Figure]^ our data read warmer 
air at 5m than at 8m. Therefore one should control the ventilation system to 
avoid energy dissipation, and/or use waste energy techniques [TSl[TB] . 

The energy balance calculations in Section 3 give the thermal radiation from 
the ceiling as the largest contribution (74%), followed by lighting (14%). To¬ 
gether with the data measurements discussed in Section such quantitative 
knowledge can aid energy saving efforts in a wide range of energy saving studies, 
since these processes occur in an average-sized ice hall under standard operating 
conditions. 

Furthermore, the analytical formula Eq.(54| holds for any problem pertain¬ 
ing one-dimensional heat conduction inside a generic sample with thermal dif- 
fusivity a, under time-dependent boundary conditions at the two surface^ 

In the bottom-up approach adopted here, all the physics is encoded in the 
boundary conditions, circumventing an otherwise involved phenomenological 
analysis. This article contains therefore a methodology which is firmly grounded 
on experimental data, using induction to obtain theoretical (predictive) results, 
and deduction to apply these formulas and to check the data consistency. 

The accurate measurements, energy balance analysis and the analytical tem¬ 
perature profile in the ice pad presented in this work can constitute useful tools 
for a proper scientific assessment, in view of today’s pressing energetic efficiency 
demands. 
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Appendix A. Analytical temperature profile during resurfacing. 

The temperature profile (54) was obtained in Section]^ in implicit form. In 
this section we put it into context, by using the initial and boundary conditions 
Eqs.(38), (40) and (41) to obtain the temperature profile inside the ice pad 
during the resurfacing process. First, we expand Eq.(51), 

f sin mr 


Sn{t) = — 
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mr 

1 

mr 


\ mr 


- 1 Ts{t) - 


sm mr 
mr 


— cos mr ) T/ (t) 


8.4- 
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then the coefficients (53) 
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because 0.998 ~ 1. Already at this stage we notice the suppression factor 1/n. 
The integral in Eq. (54) is now computed as 
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since *. 
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Figure A.16: Contributi on of the first 30 Figure A.17: Largest (n = 1) contribu- 
terms in the summation ||A.4[|, at t = 30s. tion, where 0 < t < 60s. 


We can thus recast the overall summation as 

( f sin(A„a;) 

n=l '^•^0 J 

°° f -0.01327n^ 

=eU"W+—^—10-73: 


n—l 


^ sm UTT 


— 0.01 cos riTT — 2.28 ) > sin —x 

(A.4) 


which is clearly convergent, since everything is proportional to oc l/n“. 

The first 30 terms in the summation at t = 30s are plotted separately in 
Fig A.16 we see that only the first six or seven matter significantly. The largest 
contribution, for n = 1, is shown in Fig lA.lTl where the correction to the 
temperature (in absolute value) is maximal at t = Os and minimal at f ^ 50s. 

The temperature profile inside the ice pad Eq.(54) is accordingly rewritten 
as follows, 


u(x, t) = T,ce(x, t) = [Ti{t) - Ts{t)] + Ts(t) 

+ EI / sin(A„a;) 

n=l '^•^0 J 

= -(1.4 X 10"®t^ - 0.0012t2 - O.OlOSf + 4.1438)- 

L 

+2 X 10"’^f^ - 0.0015^2 - 0.0065t - 1.0633 
+ ^ I/" sin(A„a:), (A.5) 

n=i J 


where the temperature correction generated by the summation is computed with 
dot and COI ; it is shown in Figure [A. 18| after 30s, where the first 100 terms 
are summed. 

Figure [^compares the analytical solution (A.51 to an FEM calculation for 
t = 30s, with n = 1...100. We notice a good agreement between the two curves. 

We remark that this result is specifically valid for any t < 60s, since the 
boundary conditions Ts{t) and Tj{t) are interpolations corresponding only to 
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Figure A.18: Temperature correction at t = 30s with n = 1...100 in Eq.(|A.4||. 


such time interval. Choosing different times will change the explicit form of 


both Tsit) and T/(<). On the contrary, the temperature profile Eq.(54| found 
in Section|^is given with implicit initial and boundary conditions, which makes 
it general and applicable to a range of diverse engineering problems. 
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